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Abstract 

Heavy-ion collisions provide the only laboratory tests of relativistic 
quantum field theory at finite temperature. Understanding these is a 
necessary step in understanding the origins of our universe. These lectures 
introduce the subject to experimental particle physicists, in the hope that 
they will be useful to others as well. The phase diagram of QCD is briefly 
touched upon. Kinematic variables which arise in the collisions of heavy- 
ions beyond those in the collisions of protons or electrons are introduced. 
Finally, a few of the signals studied in heavy-ion collisions, and the kind 
of physics questions which they open up are discussed. 


1 Why study heavy-ion collisions 

The universe started hot and small, and cooled as it expanded. Today vast parts 
of the universe are free of particles, except for photons with energy of about 3 
Kelvin or lower. This energy scale is so far below the mass scale of any other 
particle that no scattering processes occur in this heat bath of photons. So we 
may consider them to be free. 

This was not always so. Earlier in the history of the universe the temper¬ 
ature T was comparable to, or larger than, the mass scales of many particles. 
As a result particle production and transmutations were common. In those cir¬ 
cumstances would it be correct and useful to treat this fluid as an ideal gas? 
Such a gas cannot give rise to freeze out, phase transitions or rapid crossovers, 
and transport. We see the signatures of several such phenomena today, so we 
know that the ideal gas treatment would not work at all times. 

In the early universe many of the component particles of the fluid were 
relativistic. Since we wish to describe particle production processes in this fluid, 
we are forced to use quantum field theory at finite temperature to describe the 
contents of the early universe. The main theoretical tools required to study 


1 


thermal quantum field theory (TQFT) are effective field theory (which includes 
hydrodynamics and transport theory) and lattice field theory. Perturbation 
theory plays a limited but very important role, due to our detailed understanding 
of the technique. In order to test the formulation of TQFT, we need to think 
of experiments which can be performed easily. 

Experimental tests of TQFT in the electro-weak sector turn out to be un¬ 
feasible. Initial states made of leptons may achieve energy densities of the order 
of 1/fm 4 . However, mean-free paths due to electro-weak interactions are of the 
order 100 fm. So it is very hard to thermalize this energy density. Initial states 
of hadrons, on the other hand, have mean-free paths of the order of 1 fm, so 
the initial energy may be converted into thermal energy. By using heavy-ions, 
one can increase the initial volume significantly, and so improve the chances of 
producing thermalized matter. This is why heavy-ion collisions (HICs) are used 
to test TQFT. 

The objects of experimental study should be as many as possible, in order to 
subject TQFT to as many tests as can be conceived. The most important phe¬ 
nomena are transport properties: the electrical conductivity (important for the 
freezeout of photons), viscosity (responsible for entropy production), the speed 
of sound, the equation of state and so on. But perhaps the most interesting 
objects of experimental study are the possible phase transitions and crossovers 
associated with the symmetries of the standard model. Corresponding to every 
global symmetry there is a chemical potential. So the phase diagram of the stan¬ 
dard model has high dimensionality and potentially many phases. Experiments 
which are feasible in colliders can reach only a small fraction of the phases. 

2 Symmetries and states of QCD 

Phase diagrams display the conditions under which global symmetries are bro¬ 
ken or restored. Heavy-ion collisions explore the phase diagram of QCD. The 
global symmetries of this theory are chiral SUl (IV/) x SUr(IV/) x Ub( 1), where 
Nf is the number of flavours of light quarks, the subscripts L and R stand for 
left and right chirality, and B for baryon number. 

Chiral symmetry is explicitly broken by quark masses. QCD contains a scale, 
A QCD . Quarks with masses larger than A QCD are far from the chiral limit. The 
strange quark mass, m s , is near the scale of A QCD , and it is a detailed question 
whether treating it as nearly chiral helps in understanding the phenomenology 
of strong interactions. The up and down quark masses are much lighter than 
Aqcd and it is useful to treat them as nearly chiral. 

The resulting SUl(2) x SU/{(2) chiral symmetry is spontaneously broken 
down to SU(2) isospin symmetry in the vacuum. Signals of this symmetry 
breaking are the fact that the QCD vacuum contains a non-vanishing chiral 
condensate, (ipi/)), and that pions are massless. Departures from chirality are 
important and treated in chiral perturbation theory [T]: the most important 
result is that pions get a mass proportional to the square root of the quark 
mass. 
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Figure 1: At a critical point order parameters change abruptly, and specific 
heats and susceptibilities may have singularities. The location of the singularity 
is unique. At a crossover there are no singularities. The order parameters 
may have a large continuous change. It is possible that specific heats and 
susceptibilities peak as the temperature changes. The locations of maximum 
slope of the order parameter, or the peaks of susceptibilities, generally depend 
on the observable chosen. 


As the temperature of the vacuum is raised, keeping the baryon number 
and charge densities at zero, the condensate changes to a very small value, 
proportional to the quark mass. From thermodynamic arguments, models, and 
lattice QCD computations it is known that the change is gradual (see Figure 
[I]). One may try to characterize a temperature where this crossover happens, 
but it is a conventional number [2]. The crossover temperature, T c , depends 
upon which physical quantity is examined, but it is perfectly well-defined after 
a choice is mad«Q We make the choice that T c is given by the peak of the 
Polyakov loop susceptibility. 

The quark number for each of the Nf flavours is conserved. For the study 
of the phase diagram we need to keep in mind the up and down quark numbers 
(or, equivalently, the baryon number and the net isospin). A grand-canonical 
ensemble for QCD would then need two chemical potentials, p u and pd (or ps 
and pi), and the temperature T: so the phase diagram is three dimensional. 
As a first approximation one treats the up and down quark masses to be equal, 
and examines the phase diagram in the two dimensional slice with T and ps 
[5], and independently, of that in T and pi [4]. There has been little study of 
the more complete (and complicated) phase diagram [5j. 

The phase diagram in T and ps for small p b and non-zero light quark mass 
(see Figured]) was first investigated in [3]. At small ps the states of QCD are 
distinguished by the value of the chiral condensate. In the low T state it is 

1 Another example of a crossover is the formation of a glass by cooling of liquid silica. The 
glass transition temperature depends on what measurement one makes on the sample of the 
glass. 
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Figure 2: The phase diagram of QCD in the T-ps plane as seen with the chiral 
condensate pSJ- A line of first order phase transitions (black line) ends in the 
QCD critical point (black dot). The fireball produced in a heavy-ion collision 
lives in a track within the shaded domain. The lower edge of this domain is 
called the freezeout curve. The domain is traced out by tracks of the history 
of the fireball as the collider energy changes. For small energies, the domain 
ends at T = 0 and a chemical potential corresponding to nuclear matter. As the 
energy increases, the domain moves to small ps- The logic of the beam energy 
scan is that this domain is likely to include the QCD critical point. 


large, but becomes small at T > T c . At sufficiently high T the dependence of 
the condensate on ps is computable, and shows a gradual variation. However, 
various arguments lead to the expectation that at low T as one changes /is there 
is a first order phase transition signalled by an abrupt change in the condensate. 
The thermodynamic Gibbs’ phase rule [6j then tells us that the phase diagram 
has a line of first order transition. As we already discussed, this cannot hit the 
Pb = 0 axis or rise to T —> oo. So it must end somewhere. The end point is a 
second order phase transition, called the QCD critical point. 

The actual location of the curve of first order phase transition curve and 
the critical end point can only be predicted by a non-perturbative computation, 
i.e., through lattice QCD simulations (Tj. However, a computation at finite 
Pb requires an extension of known techniques because of a technical problem 
known as the fermion sign problem. Many such methods have been proposed, 
and many are being explored [7j. It would be fair to say that developing such 
techniques is one of the most active areas in lattice gauge theory today. 

Till now extensive computations in QCD with varying lattice cutoffs, spatial 
volumes and quark masses has been possible using only one particular method, 
involving the Taylor expansion of the pressure in powers of p b ■ As a result the 
information available until now is fairly limited, and one would hope that the 
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future brings alternative computational schemes. The current best estimate of 
the location of the critical point is [8 

T e ~ (0.94 ± 0.01)T C and n% ~ (1.68 ± 0.06)T E . (1) 

Methods have also been developed to compute the equation of state, the bulk 
compressibility, and the speed of sound in several parts of the phase diagram 
0 . 

Two more aspects of the phase diagram of QCD are interesting, but cannot 
be described here. One is the temperature dependence of the axial anomaly. 
This has been keenly investigated in recent years 0. The other is the phase 
diagram of QCD in a strong and constant external magnetic field. This has also 
generated much work recently m- 


3 General conditions in heavy-ion collisions 

I turn now to heavy-ion collisions, which is the experimental system that can test 
the computations we discussed briefly in the previous section. In this section I 
touch upon three related questions: whether thermal matter is produced, what 
its flavour content is likely to be, and how one can control the energy content 
of this matter. 

The object of study in heavy-ion collisions is the (hopefully) thermalized 
matter in the final state. In high energy colliders matter is always formed. 
In sufficiently hard pp collisions, for example at the LHC, even soft physics 
contains enough energy to create W/Z bosons, not to speak of hadrons. The 
mere production of large amounts of hadronic matter is not of interest. What 
we need to know is whether this matter re-interacts with sufficient strength to 
thermalize. In the language of particle physics this is about final state effects. 

In order to understand the time scales involved it is sufficient to run through 
a simple kinetic theory argument. Let the two-body scattering cross section be 
a. Taking the number density of particles in the final state to be n, one can 
write the mean free path as 

A oc (ncr) -1 , (2) 

If the dimensionless number 1 /(A ^/n) = a/n 2 ^ 3 = 0(1) then the mean free path 
is of the same order as the mean separation between particles. In this case, 
final state collisions are numerous, and matter may come into local thermal 
equilibrium. 

When >/S~20 GeV, we know that jets are rare. As a result, we can take the 
final state particles to be hadrons, so that er ~ 40 mb. In this case n > 5/fm 3 
may be sufficient for the final state to thermalize. This number density cannot 
be reached in collisions of protons. However, heavy-ion collisions increases n by 
some power of A, so heavy-ion collisions at this energy may thermalize. At the 
LHC, n is large, so thermalization is easier. Even high multiplicity pp collisions 
may then thermalize. The thermalized system arising from these collisions is 
the fireball which is the object of study in heavy-ion collisions. 
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This treatment is sufficient for building intuition, but a quantitative analysis 
of thermalization is more complex. The rapid expansion of the fireball implies 
that simple kinetic theory does not suffice, and the theoretical framework be¬ 
comes more complex. Some relevant references are collected here m- 

The flavour content of the fireball is needed in many analyses. Again, sim¬ 
ple arguments are sufficient to gain a quick intuition about this. The flavour 
quantum numbers of the incoming hadrons are essentially contained in hard 
(valence) quarks. At large y/S, the asymptotic freedom of QCD guarantees 
that our intuition about Rutherford scattering holds, and these valence quarks 
do not undergo large angle scattering. As a result, the incoming quantum num¬ 
bers are mostly carried forward into the fragmentation region. In terms of the 
pseudo-rapidity 



(3) 


(where 9 is the scattering angle) the fragmentation region is the region of large 
|r?|, and is called so because (classically) one finds the unscattered fragments of 
the initial particles here. 

Although the valence partons individually contain large momenta, there are 
only three of them in a baryon. Soft (sea) partons are much more numerous. As 
a result, quite a significant fraction of the energy is carried by all the soft partons 
together. These generally scatter by large angles and so stay in the central 
rapidity region (i.e., the region with \r]\ — 1). If this matter approximately 
thermalizes, then it makes the fireball which is the object of heavy-ion studies. 
The net-baryon and flavour content is small, the energy content increases with 
y/S. At high energies the central and fragmentation regions are expected to be 
well separated, i.e., one expects few hadrons in the intermediate region between 
them. 

At y/S ~1—10 GeV, baryon interactions cannot be analyzed in terms of 
quarks. In this regime the fireball may contain baryon and other flavour quan¬ 
tum numbers. The distinction between fireball (central) and fragmentation 
region may be weak. 

In the collision of point-like particles in quantum theory, the observables are 
the number of particles (or energy-momentum) hitting the detector at any angle 
9. The only control parameter is the center of mass energy, yfS. In collisions of 
extended objects, there is another control parameter: the impact parameter, b. 
This measures the separation between the centers (geometrical centers, centers 
of energy) of the colliding objects. However, b cannot directly be measured in 
an experiment. 

Instead we perform the following analysis. The total nucleus-nucleus cross 
section depends only on the energy, so one has 



( 4 ) 


Since cross sections are non-negative, the fractional cross section 
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Figure 3: The fraction of the total cross section can be used to define centrality 
classes. Different models of nuclear densities will map a centrality class to 
different impact parameter ranges. 


decreases monotonically as b increases from zero to infinity (see Figure [3]). As a 
result, an experimentally determined histogram of / would determine b uniquely, 
provided one knows the functional form of f{y/S , B). This is not yet computable 
from QCD, so one has to make models. 

The simplest, and oldest, model is called the Glauber model. In this, one 
assumes that the nucleus-nucleus collision is described by independent nucleon- 
nucleon collisions. The nucleons are distributed in each nucleus according to 
the density determined by low-energy electron nucleus collisions. Models which 
incorporate more phenomenology have also been developed; see [12| for more 
information. It has been realized in recent years that the lumpy distribution of 
nucleons in the initial state (see Figure 0]) cannot always be averaged over, but 
must be taken into account in these models. 





ylfm] 


Figure 4: Transverse energy profiles, 0.2 fm after the collision, in three models 
of initial states [T5j. The coarse-grained average of these distributions should 
give the nuclear density known through low-energy experiments. These rel¬ 
ativistic experiments capture the quantum fluctuations in the initial nuclear 
wave-function. 

As one sees, the connection between the impact parameter and the percent of 
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Figure 5: One of the definitions of centrality used by the ALICE collaboration 
m uses the histogram of measurements in the VZERO module. The experiment 
determines its correlation with the energy deposition in the ZDC. 


total cross section is indirect and model dependent. Also, when one realizes that 
the positions of nucleons inside the nuclei fluctuate from one event to another, 
it is clear that the notion of an impact parameter, and even the size and shape 
of a nucleus, are merely averaged quantities. For experimental purposes what is 
necessary is to classify events according to the degree of centrality. For this it is 
enough to define centrality by any measure which changes monotonically with 
6 ; for example, by multiplicity, zero degree calorimetry, etc.. Care is needed to 
relate these measures to each other through careful analysis of the data. One 
such analysis is shown in Figure [5] 


4 Hard probes 

One thinks of the LHC as an arena of hard QCD, i.e., of processes which convert 
the partons contained in protons into jets, heavy quarks, W/Z bosons, hard 7 , 
H and so on. The typical momentum scale in these processes is of the order 
Q ~ {x)\^S ~ 500 GeV. Final state interactions are suppressed in pp collisions 
because of two reasons. Firstly, the dense hadronic debris are separated from 
probes by large angles, A rj. Secondly, the energy scale of any remaining hadronic 
activity in the central rapidity region is small: (Et) — A QCD ~ 0.3 GeV. 

In heavy-ion collisions, the first argument can still be supported. However, 
the second argument may fail if the number density of particles, n, is large 
enough. Let us make an estimate by assuming, as before, that n = 5/frn 3 . We 
know that the actual value of n at the LHC is larger, so our argument will be 
overly conservative. Assume that the jet cone has radiutd R = 0.2, and that 
it travels about £ = 10 fin through the fireball of soft particles. Then the net 

2 The radius of a jet cone is defined to be R = Arf 2 + A <j > 2 where we take A// and A <j> to 

be the jet opening angles. 
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energy in the soft hadrons it can interact with is 

£ ~ ( E T )nRi 3 ~ 300 GeV, (5) 

where we have made a conservative estimate that the average transverse energy 
of the particles is (Et) — 0.3 GeV. Since this is comparable with the initial 
energy, final state interactions become important. An interesting consequence 
which we discuss here is jet quenching. 



Figure 6: Comparison of two-particle azimuthal distributions for central d+Au 
collisions (circles) to those seen in p+p (histogram) and Au+Au collisions 
(stars). The respective pedestals have been subtracted. 

When a jet evolves through a medium, interactions and radiation would tend 
to deplete its energy m • This simple idea is called jet quenching. The basic 
fact of jet quenching was beautifully demonstrated by the STAR collaboration 
in BNL [TB] in the plot given above. At y/S = 200 GeV jets are not very well 
developed, and one must use high-p T hadrons as proxies. STAR triggered on 
events where there is a high-p T hadron, and looked at the angular distribution 
of the next highest-p T hadron. In p+p collisions they found a peak 180 degrees 
away (see Figured]). If the trigger hadron can be assumed to come from a 
jet, then the backward peak comes from an away-side jet which balances the 
momentum. This was also seen in d+Au collisions, thus demonstrating that 
initial state parton effects in heavy nuclei do not wash away this peak. In 
Au+Au collisions they found no peak in the backward direction: implying that 
the away-side jet is hugely quencliecff]. 

A measure of the quenching is provided by a comparison of the number of 
jets of a given momentum in heavy-ion and proton collisions 

p (h ^ 1 d 3 N AA fcPN^y 1 

aa , T AA (b) dbdydp T \dydp T ) 

Here T AA is an estimate of the number of proton pairs interacting in AA col¬ 
lisions, and is usually extracted from a model, e.g., the Glauber model. The 

3 Since the near-side jet is used as a trigger, the event sample is of those in which this is 
not completely quenched. 
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numerator depends on collision centrality whereas the denominator does not. 
Energy is tremendously more likely to flow from the jet into the low-momentum 
particles in the medium (computations reveal this in phase space factors). As 
a result, one would generally expect Raa to be less than unity. 



Figure 7: The first attempt to constrain Taa from experiment. 


Since a basic input into jet-quenching is Taa-: it is important to constrain 
this through experiment. The production of high-p T photons or W/Z bosons 
provides this calibration. Since the vector bosons have no strong interactions, 
the comparison of semi-inclusive single boson production cross sections in pp 
and A A cross sections can directly measure Taa- One of the first attempts m 
to constrain this is shown in Figure [3 Small isospin corrections, shadowing, 
and initial re-scattering effects must also be taken into account more accurately 
in order to improve these constraints. 

From the observations of Raa (see Figure [8]) one can extract a measure of 
the pt change of the jet per unit distance travelled within the plasma: 



d 2 p T 

( 2^) 2 


Pt P(pt, L). 


Most attempts to extract this from data give q ~ 1-2 GeV 2 /fm, i.e., in the 
range of interest, q/T 3 ~4-5. One should be able to extract this number for 
QCD, but it turns out to be a vexing problem. There are two main methods 
to handle problems in QCD: perturbative QCD is used to compute processes 
where all momenta involved are large, and lattice QCD provides a tractable 
computational method when all momenta are small. Jet quenching couples 
a large momentum object (the jet) to low-momentum objects (the medium). 
Nevertheless there have been attempts to compute this in QCD using weak- 
coupling expansions [T9] or, more recently, lattice QCD [20]. There are also 
computations in cousins of QCD which have N = 4 supersymmetry in the limit 
of large N c [21]. 
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Figure 8: A recent compilation of data on Raa 18|. 


Raa is just the simplest of experimental variables which can be constructed. 
In order to understand how the medium steals energy and momentum from the 
jet one should also understand medium modification of rapidity and angular 
correlation, momentum imbalance between reconstructed jets, fragmentation 
functions, and jet substructure. 


5 Flow 

Once the fireball reaches local thermal equilibrium a much slower process be¬ 
gins of transport of energy, momentum, and other conserved quantities through 
the fireball. This is the hydrodynamic regime [221 . Tests of hydrodynamics 
involve the study of quantities which are called flow coefficients [23] • In order 
to understand what these are, we need to think again about the geometry and 
kinematics of the collisions. 

In the collision of point-like particles, there is a rotational symmetry around 
the beam axis. As a result, cross sections or particle production rates depend 
only on the scattering angle 0 (or equivalently, on rj or the rapidity y) and the 
transverse momentum, p T . Kinematically, there is only one initial vector in the 
center of mass (CM) frame of the problem, the initial momentum k of one of 
the particles (the other particle has momentum —k). Final state momenta see 
only the angle from k, which is 0 , and the transverse projection p T . 

In heavy-ion collisions, there is a second initial vector: b, which is the line 
between the centers of the nuclei. The existence of such a vector, not collinear 
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with k, means that the azimuthal symmetry around k is broken in the initial 
state, and final state momentum distributions may depend on angles the final 
momentum makes with both k and b as well as p T . Conventionally, these 
distributions are given in terms of 77 , p T , and the azimuthal angle (f>. The two 
vectors k and b lie on a plane which is called the reaction plane. This breaking of 
cylindrical symmetry also occurs in proton-nucleus collisions, since the proton 
can meet the nucleus with a non-vanishing impact parameter. In very high 
energy collisions, the increasing proton-proton cross section implies that the 
swollen protons can also be treated similarly. One may already be seeing such 
effects in the sample of extreme high multiplicity events in pp collisions at the 
LHC. 

The flow coefficients are the Fourier transforms of velocity distributions with 
respect to <f> [24]. The n-th Fourier coefficient is denoted by the symbol v n . 
These are normally taken at y = 0 not only because of the limited rapidity 
coverage of lreavy-ion detectors, but also because one expects the fireball to 
be well-separated from the fragmentation region. Nevertheless, studying the 
rapidity dependence of flow coefficients is of some interest. The study of the k T 
dependence of the v n is of great interest. 



Figure 9: The geometry of elliptic flow. 

Clearly, the reaction plane in different collisions can rotate around the beam 
axis, so single particle distributions will recover azimuthal symmetry when av¬ 
eraged over events. Although the overall orientation of the reaction plane is 
forgotten on the average, the relative angles between two particles remembers 
the difference from the reaction plane. So, in order to see the flow coefficients 
one has to construct the angular correlations of two or more particles. 

In the collision of symmetric nuclei, b and — b seem to be completely equiv¬ 
alent. As a result the two sides of the reaction plane seem to be completely 
symmetric. This implies that only the even flow coefficients, V 2 , V 4 , etc., are 
non-vanishing (elliptic flow is the name given to U 2 ). However, when one stud¬ 
ies the flow event-by-event (E/E) one has to take into account the fact that 
the positions of nucleons inside the nuclei may fluctuate. Then there may be 
more nucleons on one side of the plane, so breaking this orientational symmetry 
around the reaction plane, as a result of which odd harmonics may exist. Cur¬ 
rently there are studies of the directed flow v \, triangular flow V 3 , and even the 
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coefficient v$. The flow coefficients yield a combination of information on the 
initial state and the evolution of fireball. E/E fluctuations of flow coefficients 
yield more refined information on the initial state [25] 

It is claimed that the observations of v-i imply the formation of locally ther- 
malized matter in heavy-ion collisions. Although this argument is technical it 
is easy to understand this intuitively. In the off-center collisions of nuclei the 
colliding region is a pellet. Particles formed in the initial collisions have dis¬ 
tributions which have positional anisotropy, e n , the pellet being long in one 
direction (see Figure [9]). The generation of v n involves transforming e n into mo¬ 
mentum anisotropy. This is impossible unless there is hadronic re-scattering. 
Also, the measurements of V 2 show that the momentum is larger in the direction 
in which the original position distribution was squeezed. This is hydrodynamic 
flow, since that is driven by pressure gradients, and the gradients in the shorter 
direction are larger. 



Pr(OeV) 


Figure 10: Predictions for elliptic flow from ideal hydrodynamics compared to 
data. See the text for a discussion. 

The technical question here is how well does hydrodynamics explain the 
observed V 2 - Ideal hydrodynamics, i.e., hydrodynamics without dissipation, is 
a toy model which is often used to understand general features of data. This 
already come within a factor two of the data, and shows the same p T dependence 
as the observations (see Figure [TU]). It also fails in the “right” direction, in that 
it over-estimates V' 2 - Dissipation would clearly reduce the predictions, and bring 
it closer to observations [2d. One of the first results (Dusling and Teaney in 
[26] ) is shown in the figure. Intense work continues to be done in understanding 
the implications of the data and the constraints from QCD [27] . 

6 Chemical composition of the final state 

The most easily observed quantities have to do with the final state. The basic 
observables are the spectra of identified particles. The multiplicity of each 
type of particle, 7r ± , K±, etc., is called its yield. This is the integral over the 
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spectrum. Relative yields of hadrons is the outcome of hadron chemistry, i.e., 
inelastic re-scattering in the final state. Examples are, 

P + TT~ <-> Tl + 7T°, p + n~ -o- A + K°, (7) 

The rates of such predictions determine whether hadron chemistry comes to 
chemical equilibrium. 

As the fireball evolves, eventually mean-free paths or relaxation times be¬ 
come comparable to the size or expansion rate. When this happens, local ther¬ 
mal equilibrium can no longer be maintained, and hydrodynamics cannot be 
supported. Then the components of the fireball are said to freeze out. In prin¬ 
ciple freeze out could occur either before the fireball cools into hadrons or after. 
Under normal circumstances, i.e., if the thermal history of the fireball does 
not take it near a phase transition, then it seems to freeze out in the hadronic 
phase. This implies that the later stages of hydrodynamics require the equa¬ 
tion of state of a hadronic fluid. Since hadrons are massive, inelastic collisions, 
i.e.those which change the particle content of the fireball, require larger energies 
than elastic collisions. As a result, chemical freezeout, i.e., fixing of the hadron 
content of the fireball, may occur earlier than kinetic freezeout, i.e., fixing of 
the phase space distribution of hadrons. 

In the fireball the particles interact strongly enough that a temperature 
is maintained. However, at freeze out the interactions stop abruptly. So all 
hadrons emitted by the fireball at freeze out can be assumed to be an ideal gas 
of particles coming from a source whose temperature is set at the freeze out. 
This simple approximation, which goes by the name of the hadron resonance 
gas model has had remarkable phenomenological success [28]. However, recent 
measurements at the LHC (ALICE collaboration [28]) and a more careful look 
at RHIC results shows interesting discrepancies which imply that this model 
needs to be improved. 

At early times, the fireball is a reactive fluid whose description requires 
coupling of hydrodynamics with diffusion and flavour chemistry. The reaction 
rates depend on local densities as well as rates of mixing due to fluid move¬ 
ment, known as advection, as well as diffusion. In order to make quantitative 
predictions, one must first understand whether advection or diffusion is more 
important in bringing reactants together. This is controlled by Peclet’s number 

„ Lv Lv , . 

Pe = — = — = IvnM, (8) 

D t,c s 

where L is a typical macroscopic distance within the fireball over which we 
wish to compare advection and diffusion, v a typical flow velocity, £ is a typical 
density-density correlation length and c s is the speed of sound. The diffusion 
constant, D ~ £c s . We have also used the notation for the Mach number of the 
flow, M = v/c s , and the Knudsen number, Kn = L/£. When Pe C 1 diffusion 
is more rapid than advection; when Pe 1 advection is more rapid [291 . 

Peclet’s number defines a new length scale in the fireball, this is the scale at 
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which advection and diffusion become comparable 


L 


M' 


( 9 ) 


Since longitudinal flow has M < \/3, then taking £ to be approximately the 
Compton wavelength of a particle, we find that for baryons, L — 0.3 fm and for 
strange particles, L — 0.5 fm. This implies that advection may be important in 
chemical processes occurring in the early stages of the evolution of the fireball, 
but over most of its history, the availability of reactants is governed by diffusion. 

Once the reactants have been brought together we can ask whether one or 
the other reaction channel is available. If the reactions are slower than the time 
scale of transport, then we may consider the fireball to be constantly stirred. 
It is then enough to examine chemical rate equations. In this approximation, a 
toy model which takes into account only pion and nucleon reactions is: 


p = —7(p7T° — 717T + ) — ^ [pTT~ — n7T°) + • • • 

h = 7(p7r° — mr + ) + 7'(p7r - — nn°) + • • •, 

7T° = —7(p7r° — nn + ) + "/'(pTr~ — mr°) + ■ ■ ■ 

7 r+ = j(pTT° — nir + ) + • • •, 

tt~ = — 'y / (jpT^~ ~ nn°) + • • • . 


Here the label for a particle denotes the density of that particle. The rate 
constants 7 and 7 ' can be deduced from experimental measurements of cross 
sections. The equilibrium concentrations are given by 


P 

n 



7r 


(= 0 , 


( 10 ) 


where £ is the isospin fugacity. Since tt + /tt~ = ( 2 , if we set ( ~ 1, then 
pi = T log ( ~ 0. Even in this simple limit of a very rapidly stirred fireball, 
a more realistic model contains all possible reactions between many species of 
particles, of which many cross sections are unmeasured. As a result, a detailed 
model is out of reach and one must develop simplified models which catch as 
much of the physics as the state of the data justifies. 

In order to set up such a model, we consider flavour changing reactions. 
Strangeness changing processes seem to naturally split into subgroups. Indirect 
transmutations of K and n involve strange baryons in reactions such as f + 
K + -o- S° + 7T°. These have very high activation thresholds. Direct transmuta¬ 
tions can proceed through the strong interactions such as K + + K~ -o- 7r + + 7t _ . 
These are OZI violating reactions; slower than generic strong-interaction cross 
sections. Direct transmutations through weak interactions are not of relevance 
in the context of heavy-ion collisions. As a result, there is no physics forcing K 
and 7 T to freezeout together. However K and <fi are resonantly coupled, so they 
may freeze out together m- On the other hand, isospin changing processes 
(the model in eq. flOl) require extremely low activation temperatures, and may 
persist till later. 
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Figure 11: Comparison of model predictions with data in terms of the discrim¬ 
inant (model-data). The closer this is to zero, the better the model. The error 
bars show the error on the data and set the scale of what is acceptable mismatch 
between data and model. At all energies, 2CFO works better than 1CFO. 


One can capture this information into a HRG model with two freezeout 
points: one for the strange hadrons and (since this is resonantly coupled to the 
K± channel), another for non-strange hadrons. We could call this model 2CFO 
(30 in contrast with the usual HRG model with a single freeze out (1CFO). A 
comparison of measurements and best fit model predictions is shown in Figure 

El 

Interestingly, the introduction of two freeze out points allows one to do 
away with some unphysical features of the freeze out model 1CFO. In most such 
models there is a mismatch between strange and non-strange baryon production, 
which is fixed by having a fugacity factor which changes the occupancy of strange 
hadrons. This factor cannot be justified within an ideal gas picture, nor does it 
vary smoothly or monotonically as \/S is changed. Such nuisance parameters 
no longer appear within the 2CFO scheme. 

The success of the 2CFO scheme implies that as one introduces more of 
the hadron dynamics into the freezeout process, the ability to describe the data 
improves. This justifies our belief that a proper description of reactive transport 
should be able to give a good description of the final observed yields. 

The freeze out temperatures and chemical potentials in 2CFO are shown 
in Figure [T2] Also shown there is the position of the critical point of QCD 
determined in lattice studies [5]. The freeze out curves pass close to the QCD 
critical point, making it plausible that a study of the final state as one scans 
in \/~S can reveal signals of this very interesting prediction of QCD. This is 
the rationale for the RHIC Beam Energy Scan (BES) program and for planned 
future experiments in GSI and JINR. 

Experiments also measure the yields of heavy-quarkonia. In particular, the 
yield of the T family of mesons in AA collisions at LHC differs significantly from 
that in pp collisions at same \fS. This is usually reported in terms of an Raa 
for the meson. Since the quark mass is large, M>!T~ A QC d, one may expect 
that the production of quarkonia is a hard process. However the binding energy 
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Figure 12: The freeze out points obtained in 2CFO (from Chatterjee et al., 
[30]). The strange freezeout point is shown with circles, and the non-strange 
with squares. The two freeze out points at the same y/S are joined with a line. 
The large filled circle is the estimated location of the QCD critical point from 
lattice computations [8]. 


is of the order of the temperature, B ~ T, so we may expect large thermal 
effects as the cause of the change between A A and pp collisions [32] . 

Thermal lattice QCD computations show that the highest mass resonances, 
which are the least bound, are more easily disrupted at any given temperature 
[55] . This observation led to the formulation of a key observation called stepwise 
suppression, i.e., as \/S is increased Raa of the higher resonances drop below 
unity, roughly in the order of the binding energy [54] . If this works, then at a 
sufficiently high temperature it should be possible to use a thermal model to 
understand the relative yields of the T family of mesons using the variables 


r[T(n£)] 


dN r{ni) 

aiV A A 




-1 


dydpr y dydpT J 


( 11 ) 


The thermal model involves only a single parameter: the freeze out temperature 
of this family of mesons. The first data from the LHC [35] is fitted [36] well by 

Tj = 222±H MeV. (12) 

It will be interesting in future to see whether other members of the bottomonium 
family confirm this picture. Future data on r['F(n£)] will also provide useful 
tests. More detailed dynamical models m predict many more details of the 
kinematics of quarkonium suppression. 
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Figure 13: Measured (filled circles) and predicted (unfilled squares) suppression 
in the bottomonium family using a simple thermal model for freeze out [36]. 


7 Fluctuations 

Is the ensemble of heavy-ion collision events captured by a detector related to 
the ensembles required to study the thermodynamics of strong interactions? 

The least restrictive ensemble for the study of bulk matter is the micro- 
canonical ensemble. All that this requires is that the energy of a system be 
fixed. In fact, since the fireball is well-separated from the spectators, one may 
expect a mapping between the collisions and microcanonical ensemble to be 
good. However there are two obstructions, neither of them absolute. The first 
is that the microcanonical ensemble requires the energy of each member of the 
ensemble to be the same. This is hard to ensure without more control on cen¬ 
trality fluctuations than is possible at present. Secondly, one needs 47r-detectors 
to capture the entire energy of the fireball. Most detectors in use today miss a 
very large fraction of the energy. 

So one must to try to map the ensemble of events recorded in the detector 
to either the canonical or a grand-canonical ensemble. The difference between 
these is that the system being studied must exchange either energy or material 
and energy (respectively) with a much larger system called a heat-bath. Since 
detectors accept particles from only a part of each fireball, one may be able 
to map the events on to a grand-canonical ensemble. Of course, thermal and 
chemical equilibrium is necessary in order to be able to do this. 

We have already discussed the evidence that there is a degree of thermal 
and chemical equilibration at freeze out. So, if one observes a small part of 
the fireball, it may be possible to treat it in a grand-canonical ensemble where 
the rest of the fireball acts as the heat-bath. In order to make sure that the 
system (observed fraction of the fireball) is much smaller than the heat-bath (the 
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Figure 14: Selecting a volume in the fireball by selecting detector cuts works 
best when the momentum of particles in the volume are not totally arbitrary. 
Since hydrodynamics works, we know that in a small volume these momenta 
are aligned with the local fluid momentum, being smeared only by an amount 
of order T. 


unobserved fraction), one should use as small an angular coverage as possible 
while keeping the observed volume much larger than any intrinsic correlation 
volumes in the fireball (see Figure [Til) . If the acceptance in rapidity is Ay, V s 
is the size of the system, and 1 7 b is the volume of the heat-bath, then 


V s _ Ay 
V s + V b ~ 2 log (VS/M P ) ' 


(13) 


Taking Ay = 2, one sees that V b /V s is about 4.3 at the top RHIC energy of 
200 GeV, and around 7.5 at y/S = 5 TeV. These may be acceptable numbers. 
However, at yfS = 20 GeV the ratio drops to 2, and by y/S = 5 GeV, the “heat- 
bath” is smaller than the “system”. In order to keep the ratio V b /V s fixed, one 
has to decrease Ay with the beam energy. 

This may give rise to another problem, which is to keep the observed vol¬ 
ume much larger than correlation lengths. If freezeout occurs at time tj , then 
the acceptance region, Ay, corresponds approximately to a distance As = 
Tf sinh(Ay). As long as correlation lengths are linear in the inverse freezeout 
temperature 1 /Tf, it is interesting to examine 


A xTf = ( TfTf ) sinh 


( 21ogv Ss \ 

l^l + v b /v s J ' 


(14) 


If one wants V b /V s ~ 4 at yfS = 5, where Tf ~ 145 MeV, and one takes t/ ~ 5 
fm, then one finds A xTf — 2.5. This is a reasonable number, but it implies 
that Ay = 0.65 at this energy. Such a small acceptance window may cause 
statistics to drop significantly. However, for yfS > 20 GeV, there is a good 
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Figure 15: A comparison of data on fluctuations of net proton number in the 
STAR experiment with the lattice QCD computations of m from i nr . 


possibility that all these constraints may be satisfied simultaneously. Of course, 
if correlation lengths become very large at some y/S then all these arguments 
fail, and the system cannot be treated as being in equilibrium. 

Conserved quantities, such as the net particle number or energy, can change 
by transport across the boundary of the system. As a result, energy and net 
particle numbers fluctuate in grand-canonical thermodynamics. These fluctua¬ 
tions can now be mapped into E/E fluctuations. They were first discussed and 
suggested as probes of the phase structure of QCD in [55]. The experimental 
variables which allow a direct comparison of QCD predictions with data were 
first discovered in [35], and the first lattice QCD predictions were made in [ iOl . 

The existence of fluctuations means that the baryon number or energy of 
an ensemble is not a fixed quantity, but has a probability distribution. Such a 
distribution is characterized by the cumulants, [ B n ], which are defined as the 
Taylor coefficients of the logarithm of the Laplace transform of the distribution, 


20 











P(B ), of the baryon number- 


log 


dBP(B) e~ sB 


n =1 


i-sT 
n ! 


(15) 


The cumulants are related to the Taylor coefficients of the expansion of the free 
energy 139] in terms of hb simply as 

[B n ]=V' X {n \T f ,v f B )T?- 1 , (16) 

where /is) are generalized quark number susceptibilities (x^ 1 -* is the 

baryon density) |H] . As a result, ratios of the cumulants are independent of 
the factor V s . These ratios depend on the ratios of the dimensionless quantities 
x(")(2/, IIb)TJ _4 , which can be computed in lattice QCD, as demonstrated in 
[40] • Similar ratios were also discussed in [421 . 

f 

Since Tf and [i J B was already known from the analysis of yields, the exper¬ 
imental data could be compared to the lattice computation [43]. This compar¬ 
ison is reproduced in Figure [T5j The remarkably good agreement has led to 
subsequent attempts to refine the comparison. These include following up the 
suggestions in [ 40 ] that the comparison could yield a measurement of T c given 
Tf and n B [44] or the determination of Tf and fi B given T c [45] . There has also 
been a lot of work on various corrections which may need to be applied to the 
experimental data before comparing with predictions [ 46 ] . There is also sus¬ 
tained interest in fluid dynamical effects [Ml|47|. In the meanwhile, much new 
experimental data has been added from the RHIC Beam Energy Scan (BES) 
[45] . A BES-II is expected shortly. 

The long-term goal of the study of fluctuations is to understand the evo¬ 
lution of fluctuations along the freeze-out curve traced by changing the beam 
energy [35]. At higher energies one sees preliminary agreement between lat¬ 
tice predictions and experimental observations. This agreement is expected to 
break down in the vicinity of the QCD critical point because correlation lengths 
and relaxation times grow [49] . At lower energies one might expect a return to 
roughly thermal behaviour, although extracting this is fraught with theoretical 
and experimental challenges. The BES program aims to locate and study bulk 
matter near the QCD critical point. 
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